The thermal denaturation of DNA studied with neutron scattering 
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The melting transition of deoxyribonucleic acid (DNA), whereby the strands of the double he- 
lix structure completely separate at a certain temperature, has been characterized using neutron 
scattering. A Bragg peak from B-form fibre DNA has been measured as a function of temperature, 
and its widths and integrated intensities have been interpreted using the Peyrard-Bishop-Dauxois 
(PBD) model with only one free parameter. The experiment is unique, as it gives spatial correlation 
along the molecule through the melting transition where other techniques cannot. 

PACS numbers: 87.14.gk,87.15.Zg,87.64.Bx 



Deoxyribonucleic acid (DNA) is a highly dynamic 
molecule in which the base pairs, which carry the genetic 
information, fluctuate widely. This can lead to a tempo- 
rary breaking of a "closed" pair and a local separation 
of the two strands. Local openings may be activated by 
heating. At a certain temperature, local openings of the 
double helix extend over the full molecule, resulting in 
a complete separation of the two strands. This is called 
the "melting" of DNA, and may be considered a very 
rare example of a one-dimensional (ID) structural phase 
transitioE^. DNA melting attracted attention soon af- 
ter the discovery of the double helix structure"^*^ and 
was widely studied because it shows some similarity with 
DNA unwinding in the cell. There is recent renewed in- 
terest in biology due to high resolution melting methods^. 

Despite numerous attempts, the understanding of this 
transition is still a challenge. One difficulty in making 
progress is the absence of structural information at the 
transition. This is because experimental studies of DNA 
melting, made using techniques such as UV absorbance, 
circular dichroism and calorimetric studies, do not pro- 
vide spatial information. 

Diffraction techniques can provide this information. 
Indeed, the double-helix structure of DNA was solved 
by modeling the x-ray diffraction patterns from semi- 
crystalline fibre samples, measured by Franklin et al^. 
The data were analyzed for peak positions which revealed 
the double helix structure, and also that the molecules 
could have different configurational structures^. Many 
configurations are known to exist^, however the majority 
of the work to date has been on so-called 'B-form' DNA. 

A Bragg peak contains more information than simply 
its position. Analysis of the shape and width can deter- 
mine a correlation function and its characteristic length. 
The integrated intensity gives a measure of the quantity 
of the sample that scatters coherently. These quantities 
change dramatically close to a phase transition, and scat- 
tering techniques have proved to be excellent probes^. 

We have used neutron scattering to measure the tem- 



perature dependence of a strong Bragg peak from a fibre 
sample of B-form DNA. The data have been analyzed 
to extract the widths and the integrated intensities as 
the sample was heated through the melting transition. 
The structural information was compared with calcula- 
tions using an adapted form of a mesoscopic statistical 
mechanics model for DNA, known as the Peyrard-Bishop- 
Dauxois (PBD) model^. 

While the DNA molecules in a fibre are confined, large 
scale configurational changes are possible and are driven 
by a simple change in relative humidityi^'^. DNA melting 
is observable in a fibre and, given the favourable agree- 
ment between the theory and experiment, much of the 
physics of the transition can be understood from the ap- 
plication of the PBD model. 

The fibre sample was created using the "wet spin- 
ning" methodiS, with DNA extracted from salmon testes 
(Fluka), precipitated from a 0.4 M Li-salt solution. Sam- 
ples produced in this method are semicrystallinc'' with 
the DNA molecules aligned to within 5° of the fibre 
axis^^. The water content was set by keeping the sam- 
ple in an atmosphere humidified to 75 % with deuterated 
water. This ensured a B-form configurational structure, 
proven with X-ray fibre diffraction, and significantly re- 
duced incoherent neutron scattering from protonated hy- 
drogen in the sample. 

The thermodynamic behaviour of B-form fibre DNA 
was characterized by measuring the heat capacity using 
differential scanning calorimetry (DSC) at the Institut 
Neel, Grenoble, France. The samples were hermetically 
sealed in a hastalloy sample tube. The specific heat was 
measured in a Seratam Micro DSC III calorimeter rel- 
ative to an empty reference tube. The differential heat 
flux was measured as the samples were heated from 278 
K to the maximum temperature (~ 380 K) at a rate of 
0.6 K / min. Typical data are shown in Fig. [1] The 
data are sharply peaked at 360 K for this sample, which 
represents its thermal denaturation temperature. The 
transition is not reversible. 
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FIG. 1. DifTerential calorimetry signal, Cp, from a Li-DNA 
sample similar to that used in the neutron experiments, and 
the theoretical denaturation profile calculated from Pyrococ- 
cus abyssi DNA used for the analysis where k is the fraction 
of open base pairs 



After humidification, the sample for the neutron scat- 
tering measurements was placed in a Nb envelope and 
then sealed in an Al cassette using Pb wire for the seal. 
The cassette assisted in the thermalization of the sam- 
ple and ensured that the water content for the sample 
environment stayed constant throughout the measure- 
ment. The scattering was first characterized using the 
INS spectrometer at the Institut Laue-Langevin (ILL), 
France. The main Bragg peak, found at a momentum 
transfer along the fibre axis, (5||, is seen in the reciprocal 
space map in Fig. [5^. Its lattice spacing corresponds to 
the average distance between the base pairs, a « 3.4 A. 
The peak is very broad in the momentum transfer per- 
pendicular to the fibre axis, Q±. The sample can largely 
be regarded as a quasi- ID crystal, and is treated as such. 
Some DNA features are seen off the fibre axis due to its 
double helix structure, as is a weak powder ring from 
Pb at IQI ~ 2.2 A~^. Stronger rings from Al, seen at 
IQI 2.68 A-\ were monitored throughout the experi- 
ment as a check on instrument and sample alignment. 

The neutron three axis spectrometer INS, also at ILL, 
was used to measure the main Bragg peak as a function 
of temperature. This instrument was configured with 
an incident wavelength of 1.53 A ( = 35 meV). The Q- 
resolution was defined with 40' collimation before and af- 
ter the sample, and was measured by making a reciprocal 
space map of a Bragg peak from a Si single crystal. No 
energy analysis was used, thus the static approximation 
was assumed to hold for the measurements. Temperature 
control was achieved using a liquid helium cryofurnace. 

Two scans were repeated at all temperatures. Exam- 
ples are shown in Fig. [2)d. Their trajectories were cal- 
culated for nominally elastic scattering and are shown 
in Fig. [5^. Scan 1 was along the fibre axis, through 
the centre of the Bragg peak. The Q± for scan 2 was 
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FIG. 2. a) Reciprocal space map of B-form Li-DNA. The axes 
are the momentum transfer parallel (Qy) and perpendicular 
{Q±) to the fibre axis. The strong Bragg peak is observed at 
Q|| = 27r/a, where a ~ 3.4 A is the distance between the base 
pairs along the fibre axis, b) Examples of scans, fitted using 
eq. [T] The data at 299 K represent the starting point for the 
experiment. The sample is in the melting transition at 348.8 
K. The fibre structure has collapsed at 349.1 K. 



chosen such that, when Qn = 27r/a, the direction of the 
scattered beam would be perpendicular to the fibre axis. 
This type of scan has been used to measure critical phase 
transitions in low dimensional magnets^. The tempera- 
ture steps close to the melting transition were very small 
(0.1 K) and measurements at a given temperature were 
repeated numerous times to ensure thermal equilibrium 
and reproducibility. 

The data were fitted with a Lorentzian lineshape: 
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(1) 

where Iq is the integrated intensity, Qq is the peak cen- 
tre, and F is the width which is inversely proportional to 
the correlation length along the molecule. The function 
was convoluted with the instrument resolution. A sec- 
ond Lorentzian centered at Qy « 1.5 A"-'^ was needed to 
fit the scan 2 data. The amplitudes for these two peaks 
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FIG. 3. The integrated intensities and widths of the Bragg 
peaks as a function of temperature. The points show the 
experimental results, while the lines plot the theory. For the 
widths: the solid line shows the calculation for scans 1; the 
dotted line for scans 2 assuming uncorrelated components of 
structural disorder (x = 0); the dashed line for scans 2 and 
correlated disorder with x = 0.35. 



were free parameters, however their widths were set to 
be equal in the fits. Fig. [2J5 shows examples of the fits. 

The fit results are shown in Fig. [31 The normalized 
intensities for both scans followed exactly the same curve. 
A small deviation in the data at T 345 K was due to 
the discovery of an alignment issue that was found and 
rectified. The widths were qualitatively different and are 
shown for both scans. 

The samples used for the neutron and DSC measure- 
ments were produced by the same methods and appara- 
tus with nominally the same ion concentration and hu- 
midity. While melting transitions are qualitatively the 
same, the precise temperatures involved depend strongly 
on external conditions'^. Hence, it is no surprise that 
the transition temperatures differ between neutrons and 
DSC. Subsequent analysis will discuss the transition in 
terms of a reduced temperature, 9 = T /Tc, where Tc is 
defined as where half of the base pairs are open. 

The theoretical modeling of the data required the 
derivation of the scattering profile from a finite length 
of DNA, accounting for disorder due to variations in the 
base pair sequence and to thermal fluctuations. A gener- 
alization of the profile for a ID paracrystal of finite size^^ 
results in an expression for m closed base pairs: 
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2(9_lQ||X (^("'D {''^f }) ■ "^^^ Gaussian displacements, 
CTx ||, characterize the structural disorder perpendicular 
and parallel to the molecule, and may be obtained from 



conformational analysisi^ to be {^cr|^ = 0.18 A and 

^(T^)^^^ = 0.55 A. The terms, || account for thermal 
fluctuations. Along the molecule, ^|| amounts to a Debye- 
Waller correction that can be calculated from the total 
DNA mass per base pair and from the measured sound 
spee d' 'i'^ , while approximates displacements perpen- 
dicular to the molecule axis. The parameter x character- 
izes the correlation between longitudinal and transverse 
structural disorder, with x = for no correlation. 

The PBD model^ was developed to describe the melt- 
ing behaviour of DNA. It may be used to calculate 
the probability, P (m) , of having a closed cluster of m 
consecutive base pairs at a given temperature. It has 
been recently extended to model genomic DNA with no 
free parametersiS.^ and appropriate parameters from that 
work were used in the calculations here. The parameters 
were derived using sequence information from pyrococcus 
abyssi DNA and were used to model the DSC data in Fig. 
[T]with satisfactory agreement. Natural DNA has ~ 10^ 
base pairs. While P (to) may be calculated up to ^ 10^ 
base pairs, it will scale^i exponentially with m and may 
be replaced by P (to > M) = PoC"- The expected neu- 
tron scattering can then be calculated in three regimes 
using the equation: 
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The flrst term was calculated for M = 150 and the others 
for Mq — 1000. The calculated structure factors were 
then fitted with eq. [TJ Fig. [3] shows the results, and the 
experimental data, as a function of 6. 

Theory and experiment are very well matched until 
w 1, in particular for the integrated intensities where 
the agreement is almost perfect. As shown by eq. |3l 
the intensity scales with the number of closed pairs. The 
intensities are almost constant with increasing tempera- 
ture until w 0.97, due to there being a very low fraction 
of open base pairs in this temperature range. They then 
begin to fall in a smooth manner as the number of closed 
pairs decreases. This rounded, continuous behavior of the 
transition over a finite temperature range represents the 
"genomic" limit of multistep melting proposed in model 
calculations of heterogeneous DNA'^ . 

The comparison of the widths provide spatial infor- 
mation that standard observations of DNA denaturation 
cannot. Globally, the comparison between theory and 
experiment is reasonable. For 6 < 0.97 there is a qualita- 
tive difference, with the calculated widths slowly increas- 
ing with temperature while the experimental data show 
a decrease. Subsequent experiments show that this is a 
form of annealing in the sample, with the widths staying 
constant when the temperature is decreased, and hence 
is due to the fibre structure rather than DNA denaturing. 

There are no free parameters in the calculation of the 
widths for scan 1, where Q± = 0. Theory and data match 
quantitatively with a discrepancy of < 15%, validating 
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the conformational analysigJ^. The widths hardly vary 
with temperature until — 1, showing that the base pair 
openings, which start to be significant at 6* = 0.98, do 
not cause a sharp decrease in the size of the diffracting 
clusters until denaturation has occurred. Clusters of ~ 
100 base pairs remain intact well in to the denaturation 
region and then denature as a whole. 

Scans 2, with Q± ^ 0, require the introduction of x- 
Two calculations are shown in Fig. |3l The first assumes 
no correlation (x = 0) ^i^d is unacceptably lower than 
the experimental data. A moderate value (x = 0.35) 
gives a width that matches the data for 6 « 0.97, judged 
to be representative for an "annealed" sample. It would 
be interesting to test the value of x with conformational 
analysis. 

The experimental widths of scans 2 show an upturn in 
the vicinity of the phase transition (0.99 < < 1) which 
is not accurately modeled. The theory does account for 
transverse thermal fluctuations through (■^x)' however 
this effect is small as only fluctuations from closed clus- 
ters can contribute. More likely, the upturn is due to 
thermally induced misalignment of the molecules with 
respect to the fibre axis. A misalignment of 10% would 
change the width along the fibre axis by SQ\\/Q\\ = 2%, 
but the projection at non-zero Q± would create a change 
of 17%. The assumption of a ID structure begins to 
break down at these temperatures and the theory would 
have to be extended to account for that. 

The theory and experiment cease to match at 9 1 
when the experimental intensities fall precipitously while 
the theory continue to show a narrow but smooth decay. 
The data in Fig. show examples of the scattering 
below and above this temperature. The fall is due to 
the collapse of the fibre structure, as verified by optical 
microscopy on similar samples and by visual inspection of 
this sample after the experiment. The structural change 
in the fibre is distinct from DNA melting, which refers to 
the denaturing of the double helix, and is not explicitly 



described by the theory. Nevertheless, the theory gives 
some insight to the phenomenon. As is approached 
the open regions will quickly grow in size. Single strands 
are very flexible and can gain entropy by losing their 
initial orientation. The resulting structure would have 
the remaining closed, relatively rigid regions embedded 
in a liquid-like medium. Once the open regions reach a 
certain average size, the closed regions will have flexibility 
to move. The highly oriented fibre structure will break 
down, heralded by the upturn in the scan 2 widths and 
the disappearance of the Bragg peak seen in Fig. [2b- 

Despite the loss of the Bragg peak, it is interesting to 
note that the theory predicts that the size of the closed 
segments is still large (~ 80 — 90 base pairs) slightly 
above 6 = 1. This is indirectly confirmed by agarose gel 
electrophoresis analysis before and after the experiment. 
Before heating, the DNA had > 20000 base pairs, while 
after it had broken in to fragments of ~ 100 base pairs. 
The breaking can be due to the high stress concentration 
that occurs at the end of the rigid closed segments, linked 
to each other by flexible single strands, as they rotate in 
the breakdown of the fibre structure. 

In conclusion, neutron scattering has been used to give 
spatial information on the melting of B-form DNA in fi- 
bre form. The data could be understood using simple 
non-linear lattice dynamics theory with minimal free pa- 
rameters, showing the transition to be smooth and con- 
tinuous. The comparison was excellent to a point mid- 
way through the melting transition where the analysis 
was complicated by the collapse of the fibre structure. 
Results showed that large clusters of closed base pairs 
continue to exist well in to the denaturation regime, and 
showed correlations between thermal fluctuations along 
and perpendicular to the molecule. 
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